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^ ■ Abstract 

O ' 

. We use the considered axial deformed relativistic mean field theory to perform system- 

atical calculations for Z = 112 and 104 isotopic chains with force parameters NL3, NL-SH 
■ and NL-Z2 sets. Three deformed chains (oblate, moderate prolate and super-deformed 

chain) are found for Z = 112 and 104 isotopic chains. It is found that there is a chain of 
super-deformed nuclei which can increase the stability of superheavy nuclei in the Z = 112 
isotopic chain. Shape coexistence is found for Z = 112, 104 isotopic chain and the posi- 
tion is defined. For moderate prolate deformed chains of Z = 112 and 104, there is shell 
closure at N = 184 for moderate prolate deformed chain. For oblate deformed chain of 
Z = 112, the shell closure appears around at N = 176. For super-deformed chains of 
Z = 112 and 104, the position of shell closure have strong parameter dependence. There 
is shell anomalism for oblate or superdeformed nuclei. 

Key words: relativistic mean field, superheavy element, deformed configuration, shape 
coexistence, superdeformed chain 
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1 Introduction 



Since the possible existence of superheavy elements was predicted in 1960s by nuclear theo- 
reticians, the search of superheavy elements in nature has become a hot topic for scientists. 
Empirically, many heavy elements were identified by nuclear synthesis. At first, the elements 
Z = 105 ~ 108 were successfully produced and both physicists and chemists agree on the 
existence of these elements although their half-lives are not very long. During 1995 ~ 1996, 
because of the participating of more and more large laboratories in researches of new elements, 
the elements Z — 110 ~ 112 were produced by Hofmann et al. at GSI in Germany' 1-3 ' . Z = 114 
was produced by Oganessian et al. at Dubna in Russia in 1999I 4,5 '. One year later it was again 
reported that Z = 116 was synthesized at Dubna' 6 '. Recently, 288 115 and 287 115 were synthe- 
sized at FLNR, JINR^ . The achievement in producing new elements speeds up the researches 
on superheavy nuclei not only in experiment but also in theory. Theoretically, several models 
have been used and many aspects such as the collisions, structure, and stability have been 
investigated I 8-18 ' for heavy and superheavy elements. 
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t E-mail: zhongxianhui@mail. nankai . edu . cn 
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§E-mail: ningpz@nankai.edu.cn 
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Recently, there are systematically calculations for superheavy elements within the frame- 
work of deformed relativistic mean- field (RMF) theory [ 21 ~ 23 1. These calculations predicted 
that there are shape coexistence and super-deformation in the ground state of superheavy nu- 
clei and deformation can be an important cause for the stability of superheavy nuclei based on 
a constraint RMF calculation. The conclusion changes the usual conception that the existence 
of superheavy elements is due to their spherical shell structure and makes people recognize 
that deformed configurations are as important as the spherical one for stability of superheavy 
elements in theory. 

However, the constraint RMF calculation consumes much time, thus the calculation is very 
limited. Although there are many systematical calculations for superheavy nuclei which find 
there are a moderate prolate solution, an oblate solution, and a superdeformed prolate solution 
for the same element' 22 ' 23 ', there is no calculation which can give out deformed configurations 
of a whole isotopic chain. No information tells us that whether there are two or more deformed 
configurations for most of the nuclei in a whole isotopic chain or only exist in a few nuclei in 
the isotopic chain. Our aim is to study these properties in a long isotopic chain. Thus, we 
select Z = 104 and Z = 112 isotopic chains and carry out systematical calculation within the 
framework of deformed relativistic mean-field (RMF) theory. The binding energy, quadrupole 
deformation, root mean square radii, shape coexistence and shell closure are the investigative 
projects. 

2 The formalism of the relativistic mean-field theory 

The relativistic mean-field theory has been widely used to describe finite nuclei, in the RMF 
method, the local Lagrangian density is given as' 24 ' 25 ' 

C = ^{i^d^ - M)ip - g a ipa^ - gjrfu^ - g P ^P a ^ a ip 
+ \^ad,a -\m^ -\g^ -\g^ 

J-F^F^-e^YA^il-r 3 )^. (1) 

The meson fields included are the isoscalar a meson, the isoscalar-vector u meson and the 
isovector-vector p meson. M, m a , and m p are the nucleon-, the a-, the u- and the p- 
meson masses, respectively, while g ai g w , g p and e 2 /47r = 1/137 are the corresponding coupling 
constants for the mesons and the photon. The isospin Pauli matrices are written as r a , r 3 being 
the third component of r a . The field tensors of the vector mesons and of the electromagnetic 
fields take the following form: 

p^u = q^a u -O u A^. (2) 

The variational principle gives the equations of motion. For the static case the meson fields 
and photon field operators are assumed to be classical fields and they are time independent. 
They are replaced by their expectation values. The symmetries of the system simplify the 
calculations considerably. In all the systems considered in this work, there exists time reversal 
symmetry, so there are no currents in the nucleus and therefore the spatial vector components 
of uj^, p afl and vanish. This leaves only the time-like components, u°, p a0 and A . Charge 
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conservation guarantees that only the 3-component of the isovector p 00 survives. Finally we 
have the following Dirac equation for the nucleon: 

{-iaV + V(r) + (3[M + S(r)}}^ = e^, (3) 

where V(r) is the vector potential 

V(r) = g^\r) + g p r 3 p 00 + e I±^°(r), (4) 
and S(r) is the scalar potential 

S(r) = g a a(r). (5) 

The Klein-Gordon equations for the mesons and the electromagnetic fields with the densities 
as sources are 

{-A + ml}a(r) = - g a p s (r) - g 2 a 2 (r) - g 3 a 3 (r), (6) 
{-A + ml}u (r) = g„p v (r), (7) 
{-A + mJ}poo(r) = g p p 3 (r), (8) 
—AA°(r) = ep c (r). (9) 

The corresponding densities are 

A 

Ps(r) = E^W^W, (10) 
i=i 

Pv(r) = D'W^W, (11) 
i=i 

Ps(r) = E^J(r)r 3 ^(r), (12) 
i=i 

Pc(r) = E^JW((l-r 3 )/2)^(r). (13) 
i=i 

Now we have a set of coupled equations for mesons and nucleons and they will be solved 
consistently by iterations. 



3 Numerical calculation and analysis 

The validity of deformed relativistic mean-field (RMF) theory in the calculation for superheavy 
nuclei is tested in previous papers' 21-23 ' , so we do not test the validity any more in our work. In 
the process calculation, three typical sets of force parameters NL3' 19 ', NL-SHt 20 ', and NL-Z2' 10 ' 
in RMF model are chosen. The method of harmonic basis expansions is used in solving the 
coupled RMF equations. The number of bases is chosen as Nf = 12, JV& = 20. Pairing has been 
included using the BCS formalism. In the BCS calculations we have used constant pairing 
gaps A n = A p = W.lfsfA MeV' 30 l This input of pairing gaps is used in nuclear physics 
for many years. Although the BCS model may fail for light neutron-rich nuclei, the nuclei 
studied here are not light neutron-rich nuclei and the RMF results with BCS treatment should 
be reliable. The different inputs of /3q lead to different iteration numbers of the self-consistent 
calculation and different computational time, but physical quantities such as the binding energy 
and the deformation do not change much, which is tested in Ref.' 22 l Thus, when we carry out 
calculation, we only choose a proper initial fl and neglect its effect on our results. 
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3.1 Quadrupole deformation 

The quadrupole deformation of isotopic chain Z — 112, 160 < N < 200 with the force pa- 
rameters NL3, NL-SH and NL-Z2 are listed in figure 1. There are three deformed chains 1, 2 
and 3, denoted with circle, triangle and star respectively, which can be seen for all the force 
parameters NL3, NL-SH and NL-Z2 in figure 1. 1 is a oblate chain, the quantities of quadrupole 
deformation (3 2 < —0.3. 2 is a moderate or light prolate deformed chain except that there is 
light oblate deformation in the NL-Z2 calculation when N > 186. It is very interesting that 
3 is a super-deformed chain with (5 2 > 0.4 except /3 2 =0.34, 0.33, 0.32 for JV=196, 198, 200 in 
the NL-SH calculation, but they are still around /32=0.4. The phenomenon of several deformed 
configurations for a single superheavy element is predicted in RefsJ 22 ' 23 ] with constraint RMF 
calculation. In our work, it is the first time to obtain several deformed chains in a isotopic 
chain without using any constraint RMF calculations at all. 

Is it a general phenomenon that there are several deformed chains for a isotopic chain of 
superheavy elements? To answer this question, we also perform RMF calculation for Z = 104, 
152 < iV < 198 isotopic chain. The results are shown in figure 1 as well. All the results of 
the force parameters (NL3, NL-SH, and NL-Z2) show three deformed chains (oblate, moderate 
prolate and super-deformed prolate deformation ) for the Z = 104 isotopic chain as well as 
Z = 112. Although there are three deformed bands for the Z = 104 isotopic chain, the super- 
deformed configurations appear with neutron number N > 168 (for NL3 and NL-Z2 calculation) 
or N > 170 (for NL-SH calculation). The oblate deformation appears in the region of iV < 168 
for NL3 calculation, N < 170 for NL-SH calculation, however, it appears in the whole region 
of Z = 104 isotopic chain for the NL-Z2 calculation. Why oblate deformation does not appear 
in the whole region for NL3 and NL-SH calculations? Oblate deformed configurations which 
are predicted by NL-Z2 set but not by NL3 and NL-SH sets in the region. Are they physical 
solutions? In figure 1, we see there is a sudden change at N = 176 in the oblate chain for the 
NL-Z2 calculation, when N > 178 the oblate deformation changes gradually with the variation 
of the neutron number. It indicates that the results of NL-Z2 calculation in the region N > 178 
are physical solutions. The NL-Z2 set may be better than the other parameters in describing 
the deformation in larger neutron region. From the above mentioned analysis, we conclude that 
it is a general phenomenon that there are several deformed configurations for a isotopic chain of 
superheavy elements, and that the super-deformed nuclei can exist in the superheavy elements. 
Since there are several deformed configurations for the superheavy elements, is there shape 
coexistence in a superheavy element? This phenomenon will be discussed in the subsection 3.3 
in detail. 

From figure 1, we can see there are minima or maxima in the deformed chains. In super- 
deformed chain of Z = 112, the minimum appears at N = 172 for NL3, N = 168 for NL-SH and 
iV = 176 for NL-Z2, and the maximum appears at iV = 182 for all the force parameters. And 
in super-deformed chain of Z = 104, there is no minima for NL3; for NL-SH, the minimum is 
at iV = 188; for NL-Z2, sudden changes occur at N — 182, 184, which is the very minimum. In 
the moderate deformed chain of both Z — 112 and Z = 104, the minimum appears at around 
N = 184 for all the three parameters. In oblate chain of Z — 112, there is a maximum at 
N = 176 for all the force parameters. There is no obvious kink for NL3 and NL-SH calculation 
in oblate chain of Z = 104, but there is a large peak at N = 176 and a small peak at N = 190 
for NL-Z2 calculation. Generally the minimum of the prolate deformed chain and the maximum 
or peak of oblate deformed chain correspond to shell closure. In the moderate deformed chain, 
the minimum (/3 2 ~ 0) appears at N ~ 184, which agrees with the prediction in re f s J 10 > 31 ~ 35 ] 
that spherical neutron shell closures occur at N = 184. For superdeformed heavy nuclei the 
shell closures depend on the force parameters. There is strong trend that shell closure occurs 
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at iV = 176 for Z — 112 in the oblate deformed chain, for all the parameters giving a peak 
at N = 176. If the minimum or maximum is the sign of shell closure, the above analysis 
indicates that the positions of shell closure for the super-deformed and oblate deformed nuclei 
are different from that of spherical nuclei and have strong parameter dependence. We believe 
that there may be shell anomalism in oblate and super-deformed prolate superheavy nuclei. In 
fact, the phenomenon of shell anomalism is predicted in some Refs' 36 ~ 38 l 

3.2 Root mean square radii (rms) 

The root mean square (rms) radii are the project of investigation, for they contain a lot of 
important information of ground state properties. In figure 2, the rms radii of oblate, moderate 
prolate and super-deformed prolate configurations with NL3, NL-SH and NL-Z2 for the isotopic 
chains of Z = 112, 104 are listed. The solid symbol stands for neutron and empty symbol stands 
for proton radii. The circles, up triangles and down triangles denote the oblate, moderate 
prolate and super- deformed configuration respectively. 

In figure 2, it is seen that there is a gap 0.15 ~ 0.38/m for Z = 112 and 0.17 ~ 0.45/m for 
Z = 104 isotopic chain between neutron rms radii and proton rms radii of the same deformed 
configuration for all the force parameters. The gap become larger and larger with the increasing 
of the neutron number. For the three parameters NL3, NL-SH and NL-Z2, if we select the rms 
radii of NL3 as the standard, the NL-SH calculation underestimates the rms radii about 0.05/m 
and the NL-Z2 calculation overestimates the rms radii about 0.1/m in the same deformed 
configuration for both Z = 112 and Z = 104 isotopes. We also can see three obvious chains of 
rms radii for the three deformed configurations respectively in figure 2. 

For Z = 112 (160 < N < 200) isotopic chain, NL3, NL-SH and NL-Z2 calculations show 
that the rms radii of moderate prolate deformed configuration are the smallest among the 
three deformed ones, and the rms radii of super-deformed one are the largest in the region 
N < 190, namely, nuclei with large absolute values 02 tend to have larger rms radii in this 
region, however, when N > 192 the rms radii of oblate deformation are the largest. The gap 
of rms radii between the different deformed configurations is about 0.1 ~ 0.2 fm in general. 
For moderate prolate deformed configuration, there are obvious kinks in the proton rms radii 
at iV = 184, which agrees with the results in subsection 3.1 that there is shell closure at 
N=184. For oblate deformed one, the calculation of NL3, NL-SH and NL-Z2 shows kink in the 
proton rms radii at N = 176, but the NL-Z2 calculation is not very obvious. Together with 
the prediction of subsection 3.1, we are convinced that there is shell closure at N = 176 for 
oblate deformed configuration. Finally, for super-deformed configuration, the NL3 and NL-SH 
calculations show kinks at N ~ 168, 184 in the proton rms radii, and the kink appears at 
N ~ 190 for NL-Z2, there also are obvious kinks at N ~ 186 for NL3 and NL-SH at N ~ 192 in 
proton rms radii. In summary, the characters of shell closure for super-deformed configuration 
are not very obvious in the region N < 184, and have strong parameter dependence. There 
are some sudden changes at N = 198, 200 for moderate prolate deformed configuration, at 
N = 196, 198, 200 for super-deformed one with NL-SH calculation and at iV = 198, 200 for 
oblate deformed configuration with NL-Z2 calculation. The anomalism maybe come from the 
validity of the force parameter in these regions. 

For Z = 104 (152 < N < 198) isotopic chain, figure 2 shows that the rms radii of moderate 
prolate deformed configuration are larger than those of oblate one for both NL3 and NL-SH 
calculations. The gap between them is about 0.03 ~ 0.07/m. For NL-Z2 calculation, it scarcely 
shows any gap in the rms radii between moderate prolate and oblate deformed configuration 
when 152 < N < 176 in figure 2. All the calculations of NL3, NL-SH and NL-Z2 indicate 
that the radii of super-deformed configuration are larger than that of the moderate. The radii 
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of oblate deformed configuration for NL-Z2 in the region N > 178 lie between the radii of 
super-deformed configuration and the moderate one. When N > 186, it is very close to the 
super-deformed configuration's. As a whole, nuclei with large absolute values fa tend to have 
larger rms radii. In figure 2, obvious kinks can be seen at iV = 184 in the proton rms radii 
for moderate deformed configuration of all the parameters, which consists with the prediction 
in subsection 3.1 that iV = 184 is a magic number and shell closure exists there. Then we see 
the results for superdeformed configuration in figure 2. The rms radii suddenly become small 
at iV = 172 with NL3 calculation. N = 172 is a magic number for spherical nuclei in many 
Refs.[ 8 ' 10 l Is N — 172 still a magic number for superdeformed configuration, and does the 
sudden change come from the shell closure at iV = 172 as well? It is very difficult to answer 
the question, because the other two force parameters NL-SH and NL-Z2 can not reproduce 
that phenomenon. There is also sudden change at iV = 194 with NL3 calculation. The NL-SH 
calculation shows that the rms radii change gradually with the variation of neutron number 
and there is a kink at iV = 176, from which the radii increases much more slowly with the 
increasing of neutron. It is seen that the rms radii suddenly become small at iV = 182 and 184 
with NL-Z2 calculation. Is it caused by the shell closure at around N = 184? However, it is 
not predicted by NL3 and NL-SH calculations. Finally, let's see the results of oblate deformed 
configuration. We only list the results in the region 152 < N < 168 for NL3 and 152 < N < 170 
for NL-SH, because when we perform calculation with NL3 at N = 170 and with NL-SH at 
N = 172, the results suddenly change to moderate deformed configuration's. The calculation 
with NL-Z2 set shows there is sudden change at N=176. In summary, it is a very strange region 
in N = 168 ~ 178, there maybe exist complicate shell structure. 

From the analysis for the rms radii of Z — 112 and 104 isotopes, we find that there is strong 
parameter dependence in predicting the position of shell closure for oblate and super-deformed 
configuration. Shell anomalism maybe occur in the oblate and super-deformed configuration. 



Table 1: Binding energies for Z = 112 isotopic chains with calculation of NL3, NL-SH and NL-Z2 sets. B\, Bi 
and B 3 denote the binding energies of oblate, moderate prolate and super-deformed prolate configurations. N 
is the neutron number. 





NL3 


NL-SH 


NL-Z2 


N 


Bi(MeV) 


B 2 (McV) 


S 3 (MeV) 


Si (MeV) 


B 2 (MeV) 


S 3 (MeV) 


Si (MeV) 


S 2 (MeV) 


S 3 (MeV) 


160 


1946.55 


1958.84 


1953.41 


1948.74 


1962.90 


1955.63 


1946.00 


1955.60 


1952.48 


162 


1961.41 


1973.53 


1968.65 


1963.73 


1977.47 


1970.63 


1960.83 


1970.81 


1967.61 


164 


1975.69 


1987.53 


1982.56 


1978.47 


1991.75 


1984.57 


1975.06 


1984.51 


1981.87 


166 


1989.47 


2001.13 


1995.92 


1992.56 


2005.42 


1998.00 


1988.91 


1997.96 


1995.55 


168 


2002.88 


2014.09 


2009.13 


2005.83 


2017.59 


2011.23 


2002.49 


2011.19 


2008.79 


170 


2015.93 


2026.18 


2022.19 


2018.61 


2028.06 


2024.43 


2015.78 


2023.79 


2021.84 


172 


2028.60 


2037.76 


2035.05 


2031.03 


2040.07 


2037.43 


2028.69 


2036.09 


2034.51 


174 


2040.84 


2049.03 


2047.34 


2043.28 


2051.35 


2049.84 


2041.13 


2047.98 


2046.75 


176 


2052.75 


2059.90 


2058.66 


2054.83 


2062.30 


2060.72 


2053.22 


2058.94 


2058.62 


178 


2064.00 


2070.54 


2069.42 


2065.97 


2072.74 


2070.93 


2065.05 


2070.14 


2070.03 


180 


2074.84 


2080.39 


2079.65 


2076.35 


2082.60 


2080.87 


2076.62 


2080.34 


2080.93 


182 


2085.21 


2090.59 


2089.50 


2086.11 


2092.27 


2090.59 


2087.72 


2090.29 


2091.18 


184 


2094.79 


2100.68 


2099.03 


2095.33 


2101.54 


2099.97 


2097.96 


2100.84 


2100.90 


186 


2103.90 


2108.92 


2108.07 


2104.13 


2109.80 


2108.89 


2107.68 


2109.50 


2110.45 


188 


2112.76 


2116.72 


2116.79 


2112.56 


2118.12 


2117.82 


2117.09 


2117.75 


2119.70 


190 


2121.26 


2124.88 


2125.48 


2120.50 


2126.36 


2126.67 


2126.22 


2125.67 


2128.76 


192 


2129.42 


2132.67 


2134.31 


2128.10 


2134.23 


2135.18 


2135.07 


2133.36 


2137.73 


194 


2137.31 


2139.97 


2142.95 


2135.52 


2141.62 


2143.35 


2143.66 


2141.01 


2146.74 


196 


2144.84 


2146.65 


2151.09 


2142.76 


2148.55 


2151.00 


2151.76 


2148.87 


2155.49 


198 


2151.98 


2152.87 


2158.68 


2149.67 


2152.87 


2158.97 


2159.33 


2156.53 


2163.79 


200 


2158.78 


2165.98 


2165.99 


2156.29 


2161.52 


2166.40 


2167.00 


2163.21 


2171.58 
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Table 2: Binding energies for Z = 104 isotopic chains with calculation of NL3, NL-SH and NL-Z2 sets. Bi,B 2 
and B 3 denote the binding energies of oblate, moderate prolate and super- deformed prolate configurations. N 
is the neutron number. 





NL3 


NL-SH 


NL-Z2 


N 


r>\ ^ivie V ) 


jj2 l^-Lvie v j 


-D3 ^ivie v ) 


i)i ^ivie v ) 


jj'2 i^ivie v ) 


jr>3 ^ivie v j 


_Di ^ivie v j 


r>2 i^ivie v ) 


£>3 ^ivie V J 


152 


lo / 0.4o 


looy.zo 




1 S7fi QQ 


1 fiOO £7 

loyz.o i 




1 fiV^ OO 

lo / o.zy 


1 QQ7 nQ 

loo / .Uo 




154 


T SOD 1 A 

loyu. 14 


1 ono on 
lyuz.yu 




1 coo nc 
loyz.Uo 


lyuo.oi 




1 QQO P.A 
looy.D4 


t onn fiQ 
lyuu.Do 




156 


lyuo.zy 


1(11 C QO 

iy lo.oy 




1 ori/i 07 
iyu4.y / 


1 Ol O AQ 

iy ly.uo 




1 ono oo 
lyuz.yz 


1 Q1 Q Of^ 

lyio.yo 




158 


lyio.yo 


1 OOQ QK 

lyzo.oo 




1 oi 7 Aa 

iyi / .4d 


1 QQ 1 Q e; 

lyoi.oo 




1 Ol K CO 

lyio.ou 


1 OOft 77 




160 


i noc 1 q 


iy4u.zo 




lyzy.oo 


1 O/l Q OO 

iy4o.zz 




iyz i . ( O 


1 OQQ 71 

iyoo. I 1 




162 


i oqq on 

lyoy.yy 


mf;i /in 

iyoi.4u 




1 O/l 1 QC 

iy4i.oo 


1 flKQ OO 

lyoo.yy 




1 O/i n nn 
iy4U.UU 


i oc;n no 




164 


1 1 /I /I 

iyoi.44 


1 Ofi1 1 e; 

lyoi.io 




1 OKO QA 

iyoz.o4 


1 Ofi/I 01 

iyo4.zi 




KiKo nn 

lyoz.uu 


t ofin no 

lyou.uy 




166 


1 QftO QO 

lyoz.oz 


1 ovn 77 

iy / u. 1 1 




1 GfiQ O/l 

iyDo.y4 


1Q7/I QO 

iy / 4.0Z 




1 OfiQ 7G 

iyoo. / y 


T OfiO 7/1 

iyoy . 1 4 




168 


1 070 fin 
iy * z.ou 


1 osn 
lyoU. oZ 


1 07ft OS 


1 07/1 fifi 

iy / 4. Do 


iyo4.ui 


1 07fi A 3 
iy / 0.4o 


iy / o.4o 


1 070 no 
iy / y .uz 


iy /o.UU 


170 




1 ACQ Qfi 

lyoy.oo 


1 Ofift O/l 

iyoD.y4 


1 not; m 

lyoo.ui 


1 OOQ 1 Q 

lyyo.io 


1 OfiQ /I Q 

iyoo.4o 


1 Ofifi so 

iyoo. oz 


1 QQQ 7Q 

iyoo. / z 


lyoo.u / 


172 




1 QQQ 94 

lyyy.z^t 


1 QQ7 Qft 

lyy / .yo 




ZUUO. 14 


1 QQ7 Qft 

iyy i .yo 


1 QQ7 7R 

iyy I . lo 


i QQQ 
lyyo.oo 


1 QQ7 77 

iyy t.it 


174 




onnc f^Q 


ZUUO.Ol 




om o qo 

ZUlZ.oZ 


ZUUo.oo 


onns i /i 

ZUUo. 14 


onns no 

ZUUo.UZ 




176 




2017.68 


2014.23 




2021.20 


2015.58 


2018.38 


2017.68 


2016.55 


178 




2026.58 


2022.61 




2029.66 


2023.68 


2024.41 


2027.31 


2025.90 


180 




2034.82 


2030.82 




2037.56 


2031.36 


2033.09 


2035.74 


2034.88 


182 




2042.81 


2038.60 




2045.37 


2038.75 


2041.30 


2044.36 


2040.45 


184 




2050.54 


2045.89 




2052.65 


2045.97 


2049.34 


2052.73 


2048.21 


186 




2056.40 


2052.80 




2058.30 


2051.94 


2056.83 


2059.16 


2056.42 


188 




2061.92 


2059.48 




2063.67 


2058.63 


2063.97 


2065.98 


2063.95 


190 




2069.40 


2066.02 




2069.59 


2065.06 


2071.20 


2073.86 


2071.44 


192 




2075.42 


2072.37 




2075.05 


2071.31 


2078.14 


2080.90 


2078.79 


194 




2081.42 


2079.81 




2079.84 


2077.46 


2084.85 


2087.56 


2085.65 


196 




2087.61 


2085.75 




2086.46 


2083.17 


2091.27 


2093.73 


2092.12 


198 




2093.75 


2091.22 




2092.22 


2098.35 


2097.51 


2100.33 


2098.32 



3.3 Binding energy 

In the above subsections, we have discussed the deformation and root mean square radii of 
Z — 112, 104 isotopic chains, and find that there are three deformed chains for each isotopic 
chain. Since there are multi-deformed configurations for the superheavy nuclei, is there shape 
coexistence in them? In this subsection, we will discuss it in detail. The binding energies of the 
three deformed chains with the force parameters NL3, NL-SH and NL-Z2 are listed in table 1 
and table 2 for Z = 112, 104 respectively. B±, B 2 and £> 3 denote the binding energies of oblate, 
moderate prolate and super-deformed prolate configurations respectively. The binding energies 
with NL3, NL-SH and NL-Z2 in the same deformed chain are very close, the difference between 
them is not larger than 0.05 percent (about 10 MeV). One of the important conditions for the 
shape coexistence is that the deference of binding energies between two deformed configurations 
is very small, usually less than 1 MeV. If it is larger than 1 MeV, the translation between two 
different configurations is very difficult and the probability of shape coexistence is very little. 
To find the sign of shape coexistence, we compare the binding energies of different deformed 
chain with each other, and list the results in figure 3. The circle denotes the values of Bi — B 2 , 
triangle is for B 2 — B3 and star is for B\ — B%. In figure 3, it is seen that all the values of 
Bi — B 2 , B 2 — £> 3 and B\ — B 3 are in the region —14.5 ~ 8 MeV. It is obvious that the abstract 
values of them decrease with the increasing of the neutron number until to a certain region 
(N ~ 184) and the trend changes when iV > 184 for Z — 112 isotopic chain. For Z = 104 
isotopic chain, the phenomenon is less clear than that of Z = 112. It is interesting that these 
values are very close to the fission barriers as high as 8 ~ 12 MeV around the double-shell 
closure Z = 1U,N = 184l 27 ~ 29 l. This means the deformation is important for the stability of 
superheavy nuclei. From figure 3, we can also clearly find in which region the shape coexistence 
exists and what kind of shape coexistence is in the long isotopic chain. 
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For Z = 112 isotopic chain, from table 1 and figure 3, we can see B 2 is always the largest 
one and B\ is the smallest among B\,B 2 and B3 in the region iV < 186 for NL3, N < 188 
for NL-SH and N < 178 for NL-Z2 calculation respectively. It indicates that the moderate 
prolate deformed configuration is the ground state of nuclei N < 186 for NL3, N < 188 for 
NL-SH and iV < 178 for NL-Z2 calculation and the oblate or superdeformed configuration is 
the exciting state of them. The most distinct character is that B 3 is the largest one when 
N > 188 for NL3, N > 190 for NL-SH and N > 180 for NL-Z2 calculation, namely, there is 
a chain of super-deformed configurations which become the ground state and are more stable 
than the other deformed configurations. Although there are some differences with the different 
force parameters in predicting the region of super-deformation, all the parameters' calculation 
show there is a more stable super-deformed region in the Z = 112 isotopic chain. It agrees with 
the Bohr and Mottelson' suggestion' 26 ! that deformation can increase the stability of the heavy 
nuclei. But it should be noted that it is not for all heavy nuclei that deformation can increase 
the stability. This phenomenon appears only in special region of superdeformed chain. Figure 3 
also shows that the results with NL3, NL-SH and NL-Z2 set are different in some details, but the 
total trend is consistent. In figure 3, the star is below and far from the dotted line 8E—-1 MeV, 
it means that it is impossible for oblate and super-deformed prolate deformation to coexist in 
Z = 112 isotopes . There is no shape coexistence in the region N < 172, for no symbols in the 
region 1 MeV < SE < — 1 MeV as shown in figure 3. The NL3 results show there maybe exist 
moderate- and super-deformed shape coexistence in the region 176 < N < 190 of Z = 112 
isotopes, for the value of B 2 — B s (triangles) is around 1 MeV or within 1 MeV. Although the 
NL-SH results for B 2 — B 3 (triangles) are larger than 1 MeV in the region 174 < N < 184, they 
are very close to 1 MeV, together with the values within 1 MeV, we believe that moderate- 
and super-deformed shape coexistence maybe occur in 174 < N < 192. The NL-Z2 results 
predict the moderate- and super-deformed shape coexistence should occur in 174 < N < 186. 
It is interesting that all the regions of shape coexistence predicted by the parameters NL3, 
NL-SH and NL-Z2 contain the magic number N=184. It maybe have some trends that shape 
coexistence occur usually at around the magic number. NL3 calculation also shows moderate 
prolate and oblate deformation coexist at N = 198, Z = 112, and NL-Z2 shows moderate 
prolate and oblate deformation coexist at N = 188,190. The calculation of Ref.^ 23 ' predicts 
N = 172, Z = 112 exists shape coexistence, but our result does not indicate there is any shape 
coexistence. 

On the other hand, for Z = 104 isotopic chain, from table 2 and figure 3, we can find 
there is no super-deformed configuration until N > 168 and no oblate deformed configuration 
when N > 168 for NL3, N > 170 for NL-SH calculation. B 2 is always larger than B% except 
N = 198 for NL-SH calculation. It may be an anomalism or un-physical solution. Thus, 
for Z = 104 isotopic chain, the moderate prolate deformed configuration may be the ground 
state, and the deformation can not increase the stability in this case at all. In figure 3, the 
results of NL3 and NL-SH are very close, the triangles, circles or stars are not in the region -1 
MeV< SE < 1 MeV. It indicates there is no shape coexistence for Z = 104 isotopes with NL3 
and NL-SH calculations. But the results that NL-Z2 parameter shows are very different from 
those of NL3 and NL-SH sets. It predicts moderate prolate and oblate deformation maybe 
coexist at N = 172, 174 and 176, moderate prolate and super-deformation maybe coexist in 
168 < N < 180 and oblate and super-deformation maybe coexist at iV = 172, 174 and in the 
region 182 < N < 198. The shape coexistence also occurs around the magic number iV = 174 
and 184. 
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4 Summary 



We use the considered axial deformed relativistic mean field theory to perform systematical 
calculations for Z — 112 and 104 isotopic chains with force parameters NL3, NL-SH and NL- 
Z2 sets. First, three deformed chains (oblate, moderate prolate and super-deformed chain) 
are found for Z = 112 and 104 isotopic chains. Second, a super-deformed chain can be the 
ground states for Z = 112 isotopic chain when N > 188 with NL3, N > 190 with NL-SH and 
iV > 180 with NL-Z2 calculation. This confirms that the deformation can increase the stability 
of some superheavy nuclei. Third, although shape coexistence is a general phenomenon for 
superheavy nuclei, it only appears in some special regions of an isotopic chain, and only some 
special kinds of shape coexistence can exist. For Z — 112 isotopic chain, it is predicted that 
moderate prolate and super-deformation coexist in the region 174 < N < 184 for NL3 set, 
174 < N < 192 for NL-SH set, 174 < iV < 186 for NL-Z2 set. There is no other kind of shape 
coexistence for Z — 112 isotopic chain except moderate prolate and oblate deformation at 
N = 198 for NL3 set, moderate prolate and oblate deformation at N = 188, 190 for NL-Z2 set. 
For Z = 104 isotopic chain, there is no shape coexistence with NL3 and NL-SH calculations. 
However, the NL-Z2 calculation predicts that moderate prolate and oblate deformation maybe 
coexist at N = 172, 174 and 176, moderate prolate and super-deformation maybe coexist in 
168 < N < 180 and oblate and super-deformation maybe coexist at iV = 172, 174 and in the 
region 182 < N < 198. We predict that shape coexistence maybe occur around the magic 
number. Finally, it is found that there is shell closure at iV = 184 (f3 2 ~ 0) for moderate 
prolate deformed chains of Z = 112 and 104, at N = 176 for oblate deformed chain of Z = 112, 
the position of shell closure have strong parameter dependence for super-deformed chains of 
Z — 112 and 104. It is confirmed the prediction that there is shell anomalism for oblate or 
superdeformed nuclei. 
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Figure 1: The parameter of quadrupole deformation /3 2 for Z = 112 and 104 isotopic chains with calculation 
of NL3, NL-SH and NL-Z2 sets. Empty circle, empty triangle and empty star denote oblate, moderate prolate 
and super deformation respectively. 
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Figure 2: The root mean square radii for Z = 112 and 104 isotopic chains with calculation of NL3, NL-SH 
and NL-Z2 sets. Circle, up-triangle and down-triangle denote oblate, moderate prolate and super deformed 
configuration respectively. The solid symbol stands for neutron and empty symbol stands for proton 
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Figure 3: The energy differences for Z = 112 and 104 isotopic chains with calculation of NL3, NL-SH and 
NL-Z2 sets. Bi,B 2 and B 3 denote the binding energies of oblate, moderate prolate and super-deformed prolate 
configurations. The circle denotes the values of Bi—B 2 , triangle stands for B 2 — B 3 and star stands for Bi — B 3 . 
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